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1. Introduction 


For the topology and shape optimization of structures, the different realizations of 
homogenization method are widely used (Arora, 2004; Bendsoe &e Sigmund, 2003). This 
method is highly effective for shell constructions and allows implementing topometry and 
topography, sizing and shape as well as freeform optimization (Vanderplaats, 200) 
However, it is a very time consuming procedure because the number of design parameters 
can reach a million and more. There is the possibility of taking into account some 
technological factors, nevertheless, in the case of bulky bodies it frequently produces shapes 
that are difficult to manufacture. As shown in work (Mallerschon et al., 2010), the Hybrid 
Cellular Automata method does not allow parallelization of computations and PBS queuing 
system has been used. At the same time the following resource saving approach 
anushevskis et al, 2010; Janushevskis & Melnikovs, 2010) can be used for shape 
‘optimization which includes the following steps: 1) Planning the position of control points 
fof NURBS (see, for example, Saxena é Sahay, 2005) for obtaining a smooth shape. 2) 
Building geometrical models using CAD software in conformity with design of experiment. 
3) Calculation of responses for a complete FEM model using CAE software. 4) Building 
metamodels (surrogate models) for responses on the basis of computer experiment. 5) Using 
metamodels for shape optimization. 6) Validating the optimal design using CAE software 
for the complete FEM model, 


2. Basics of approach 


Latin Hypercube type experimental designs first proposed by Vilnis Eglajs in his work 
(Audze & Eglajs, 1977), then by McKay (McKay et al, 1979) and used by many other 
investigators, as well as its improvements (Auzins, 2004) are a very essential aspect utilized 
in the proposed method. The significance of approximations for the solution of optimization 
Problems proposed by Lucien Schmit in his early works (Schmit & Farcshi, 1971) and 
nowadays generally recognized as the Response Surface Method is also the foundation of 
current approach. The use of planned computer experiments and the metamodeling 
(surrogate model) approach (Forrester et al., 2008; Sacks etal, 1989) ensures great economy 
of computing time, especially for finite element (FEM) calculations. First of all let us 
demonstrate our approach on the simple test problems. 
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2.1 Test problem of plate bending 
A clamped square plate is considered under a concentrated load of 500 N applied at centre 
in a direction orthogonal to its main surface. The isotropic material properties are: the 
Young's modulus E = 200 GPa, the Poisson's ratio v= 03 and dimensions are 
(400x400x4.2118 mm. The shape optimization of the plate with constant thickness is carried 
out to minimize its volume in the case of a single displacement constraint 6 = 0.5 mm. The 
cutout shape of the plate is defined by subsequent techniques shown in Figure 1: 1) with the 
‘Points that are connected with straight lines; 2) with the NURBS knot points; 3) with the 
control points of NURBS polygon. Due to symmetry only W ofthe plate is considered for 
‘cutout definition and V4 of the plate for problem solution by FEM- 


D ) 


Fig- 1. Techniques for definition of cutout. 


‘Three parameters are stated to define location of points. Parameters are varied in the 
following ranges: 100s X1170; 100<X2<210; 100X3<230 mm for the first two variants and 
100¢X1180, 100X225; 100£X3<230 mm for third variant of definition. In the last case at 
‘both end points two continuity vectors are defined additionally with direction normal to the 
side and to symmetry axis of the corresponding plate and with fixed length of 19 and 3 mm. 
The design of experiment for 3 factors and 40 trial points is calculated with mean-square 
error (MSE) criterion (Auzins & Janushevskis, 2002; Auzins et al., 2006) value 0.4262 using 
EDAOpt (Auzins & Janushevskis, 2007) - software for design of experiments, approximation 
and optimization developed in the Riga Technical University. The geometrical models are 
developed using SolidWorks (SW) for all variants. The shapes for the third variant of 
definition are shown in Figure 2 In the next step responses of these models are calculated by 
SW Simulation (Lombard, 2008), using elements with a global size 4 mm and total number 
‘of DoF ~100000. Then these responses are used for approximation by EDAOpt. For example, 
for approximation of response y by quadratic polynomial the following expression (se, for 
example, Auzins & Janushevskis, 2007) is used: 


9= As SpE È pas È pae o 


Where there are d variables xy... xo L= (d+1)(d+2)/2 unknown coefficients land the errors 
‘care assumed independent with zero mean and constant variance 2. 
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Fig 2 Shapes of cutout in compliance with design of experiment. 


In the case of locally weighted polynomial approximation, coefficients fe. Ps Ph) 
depend on point xə where prediction is calculated and are obtained using the weighted least 
squares method: 


F=argmin Y wra- =H)? G] 


‘where the significance of neighboring points in the set Ns is taken into account by Gaussian 
‘kernel (weight function): 


w(u) sexp(-në) o 


Where u is Euclidian distance from x» to current point and æ is a coefficient that 
characterizes significance. 


Quality of approximation is estimated by relative crosvalidation error using leave-one-out 


crosevaldaton: 
PA -uP 
Oy, =100%. a 


Where root mean squared prediction errr stands in numerator and mean square deviation 
of response from its average value stands in denominator, n is the total number of 
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experimental tals and $f) ences approximated value for response in Hh point 
‘aula without king to acount the A, aipania pint 


Using the obtained locally weighted polynomial approximations by global search procedure 
\anushevskis et. al, 2004), implemented in EDAOpt, the optimal cutout shape is obtained 
(see Figure 3) for the different aforementioned techniques. In table 1 the results are 
summarized and compared with volume obtained in work (Liang et al, 2001) using the 
homogenization method. Variants correspond to shapes shown in figure 3, The value of 
Gaussian kernel parameter 1 of the local quadratic polynomial approximation is chosen to 
minimize relative leave-ane-out cressvalidation error ay of approximations of appropriate 
responses, ie. deflection 5 and volume v of plate. isthe predicted volume calculated using 
approximations and z» is the actual volume calculated using a geometrical model. and z 
in't show comparison of appropriate volume in respect to the volume obtained in (Liang et 
al., 2001). Best results are achieved with the technique using the control points of NURBS 
Polygon. This allows reducing the volume of the plate by 1.38 % (Fig. 3 e) in comparison 
with the homogenization method. 


i she ale a 
hog og: 


ig Shapes af plite oboe by (a) nmogeteatn meh fang 101 ely 
the current approach using different techniques: (b) the points that are connected with 
straight lines; (c) with the NURBS knot points; (d) with the contral points of NURBS 
polygon; (e) same as “d” but with additionally optimized tangent weighting at the spline 
endpoints; (f) same as “e” but with circle added. 


n GE SS um non ti x 
7 ee 5 pE z = 
bo a7 98 O0 6478 6933158 1 os 
e 17 98 O0 6898867  6S81588 O4 009% 
a 156 981 00S 6886232 6872198 016 -D0 
e 32 0 016 6777524 67800975 135 138 


‘Table 1, Quantitative indices of the shape optimization of cutout far the plate bending problem. 
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1t should be mentioned that the predicted volume for variant e is in good agreement with 
the actual value. At the same time the total number of FEM problem calculations (ie. 
number of trials of computer experiment) is less on 10 than by using the homogenization 
‘method. The distribution of resultant displacement for plate of optimal shape (variant e) is 
shown in fig 4. 


Fig 4. Displacement of plate for variant e 


2.2 Stretched plate test problem 


A square plate with dimensions 1300x13000.0001 mm is considered under two axial 
stretching loads of p = 0.65 N/mm? and p/2 accordingly (see Figure 6). The isotropic material 
Properties are: the Young's modulus E = 210 GPa, the Poisson's ratio v= 0.3. The shape 
optimization of the plate with constant thickness 3s carried out to minimize its volume in case 
of a single constraint on maximal level of equivalent stress er = 438 MPa. Due to symmetry 
anly 14 of the plate is considered for cutout definition and for problem solution by FEM. The 
cutout shape of the plate is defined by five coordinates of points 2, 3, 4, 5 and 6 situated on 
Straight lines that make angles of 0; 225; 45 and 90 degrees with horizontal axis as shown in 
Figure 5. The initial cutout shape is a circle with radins 250 mm and the volume of plate n is 
3734 maw. 


Fig.. Scheme of 4 of plate with initial cutout 
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‘The aforementioned five parameters are varied in the range from 250 to 640 mm. The design 
of experiment for 5 factors and 112 trial points is calculated with MSE criterion value 053%, 
In the next step responses of these models are calculated by SW Simulation using shell 
elements with global size 7.5 mm and total number of DoF ~10S000, Then these responses 
are used for approximation by EDAOpt Using obtained locally weighted polynomial 
approximations by global search procedure, the optimal cutout shape is obtained (see Figure 
6) for the different aforementioned techniques. In table 2 the results are summarized and 
compared with the volume obtained in work (Papadrakakis et al, 1998). Variants 
correspond to shapes shown in figure 6. 


4 
a) b) 
f y k 4 
9 a 
Fig. 6. The obtained plate shapes. Plate shapes obtained (a) in (Papadrakakis et al., 1998); and 


by current approach using different techniques: (b) the points that are connected with straight 
Jines; (c) with the NURBS knot paints; (d) with the control paints of NURBS paly gon. 


> Te E o Come Game 
a 38 % moe ma O M MP 
7 = i 7 E = ae 
b B 873 oo 2598 S a 52 
e 126 a48 27 76 D AB a3 
a 730 m% 17 00 250690 438 aM 


Table 2. Quantitative indices of the shape optimization of cutout for the plate stretching 
problem. 


The value of Gaussian kernel parameter e of the local quadratic polynomial approximation 
is chosen to minimize relative leave-one-out crossvalidation error oi of approximations of 
appropriate responses, ie. maximal equivalent stress anu and volume » of plate. u is 
predicted volume calculated using approximations and su is actual volume calculated using 
full FEM and geometrical models, Again the best results are achieved with the technique 
using the control points of NURBS polygon. This allows reducing volume of the plate by 
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105 4 (Fig, 6 d) in comparison with the volume obtained in (Papadrakakis et a, 1998). 
Distribution of equivalent von Mises stress for case d is shown in Figure 7. 


ven Maes ant) 


‘ig. 7. Equivalent von Mises stress levels in the plate for variant d. 


3. Tensometric wheel pairs 


[At the present time special tensometric wheel pairs are used for the wheel - rail system 
‘monitoring, For each type of rolling stock these wheel pairs must fit the vehicle's wheel 
‘wearing condition, diameter and bearing box connection type. Using and delivering the 
tensometric wheel pairs is expensive and takes a lot of time for preparing strength — 
dynamics tests. In this work removable equipment for monitoring is proposed for mounting 
fn the ordinary wheel pairs. The monitoring wireless system (Grigorov, 2004; Hart, 1986) for 
80 tons wagon (freight car) is taken for prototype. The movable part of the equipment (Fig. 
8) consists of a removable disk, two transmitters and a transmitting antenna as well as strain 
gauges bonded to the wheel at defined places. The removable disk is fixedly attached to the 
‘wheel pair's axis. A circular transmitting antenna and two transmitters are mounted on the 
outside of the disk. 

Removable equipment must be lightweight to minimize distortions of measurements and at 
the same time it must be with appropriate durability. During testing dynamical loads 
caused by rail joints, railroad switches and other irregularities as well as due to defects of 
Wheel geometry are transmitted from wheel pairs to the removable disk which is rigidly 
mounted on the wheelset axis. Therefore shape optimization of the mounting disk that is the 
main heavy-weight part of the equipment is very important to reduce its total weight. 
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Disk Flat of wheel 
Transmitter 


© 
Fig- 8. (a) Removable dise with elements of measurement system and (b) its mounting place. 


3.1 Loads acting on wheel pairs 


‘The main loads are acting in vertical direction and are caused by railroad irregularities and 
‘wheel defects in the wheel - railroad contact. The removable disk sustains all loads from the 
wheel pair because it is rigidly fastened. Strength of the removable disk is calculated for 
maximal possible loading. For example, in the case when the wheel pair has 2 mm flat of 
‘whee! (Fig b), the loaded and empty wagon wheelsets undergo different loads in vertical 
direction at different velocities (see Table 3) (Sladkowsky & Pogorelov, 2008). Fer an empty 
Wagon the maximal load is at velocity 5 m/s, but for the loaded wagon at 10 m/s. Strength 
of the removable disk will be analyzed with maximal vertical load - 620.6 KN for two cases 
ff orientation of the transmitters - horizontal and vertical when the wheel pair acceleration 
can reach 12 g (the gravitational acceleration g = 9.81 m/? is taken into account) 


Taima load in moment of shock KN 
Velocity of wagon m/s aoa —— 
Satie Toad ELJ 103 
T Tet Bis 
z Tz EZ | 
5 EA 372 
EI Tas Tae 
EJ EII C ] 


Table 3, Load versus velocity of wagon (Sladkowsky & Pagorelov, 2008) 


Strength of the disk under centrifugal load will also be analyzed at maximum vehicle 
velocity 200 km/h (wheel angular velocity = 11698 rad/s). In this case the disk is 
considered as new without wear on riding circle. 

Besides, frequency analysis was made to find natural frequencies of the wheel pair and 
evaluate possible resonance in the case of flat of wheel, Obtained results show that 
excitation frequencies at velocities of operating conditions are significantly smaller than 
fundamental frequency. 
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3.2 Disk model tor strength calculation 


The geometrical model of the disk is created using SW. It takes into account the shape, size and 
material of the disk and transmitters (dimensions 20x55x80 mm, mass 0.1 kg). The transmitting 
antenna (Fig. 8) is removed from the calculation model because it has small dimensions and 
lightweight material. The calculation model also doesn't consider fastening holes and 
fastening elements, The transmitting antenna works stable if its displacement in axial direction 
i less than 2 mm (Hart, 1986). This constraint is taken into account in the next optimization. 


‘The strength calculations are performed using SW Simulation. First, the shape of the 
removable disk (radius R = 300 mm and thickness b = 10 mm) is changed to an ellipse with 
semi-major axis length 300 mm and semi-minor axis E1 (Fig. 9 a) by simple size optimization 
‘of EL. This shape is convenient for the equipment mounting purposes and is taken for initial 
design. 


FE mesh (Fig. 9 b) is generated with second order tetrahedral solid elements and is 

generated to get high accuracy results. It consists of about 51000 elements with 86000 nodes 

{258000 DOF), 

Removable disk's displacement is restrained on its cylindrical face (Fig. 9 c). The material of 

the disk is aluminum alloy (1060 H12) with elastic modulus E = 69000 MPa, Poisson's ratio 
33, mass density p= 2700 kg/m’ and yield strength a, = 27.5742 MPa. 


E1190 
e 
© ® 


© 


element mesh of the model; (e) Scheme 


Fig. 9. (a) Ellipsoidal disk; (b) Computational fini 
of disk fastening. 


The material's ultimate fatigue resistance is calculated as (State Railway Research Institute 
ISRRIJ, 1998): 

1040, © 
So stresses must be less than o, =11.0297 MPa. Additionally the value of the factor of 
safety FOS = 275 is assumed to be sure that the disk will be durable in any worst-case 
situation (SRRI, 1998). The acceptable stress in the disk material is reduced to apu, =4MPa. 
Von Mises yield stress criterion is used for all strength calculations: 
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Where @;,03,2aFe principal stresses. 


‘Thereby the von Mises stress at any point of the disk should be less than acceptable stress: 
Stes “Fie o 


3.3 Stresses in initial design disk 


‘Thrve variants of stressed state of the disk are analyzed, i, from loads due to flat af wheel 
in two cases of the disk orientation: when the major axis of ellipse is vertical and horizontal 
as well as from centrifugal loads, 

We consider a loaded wagon with maximal loading in moment of shock that occurs at 
velocity 10 m/s (Table 3). Maximal stresses in moment of shock (acceleration a= 119.3 m=) 
are shown on Fig. 10. As we can see, values of maximal stress levels for both orientations of 
the disk are very similar, 


@ 
ig. 10, Von Mises stresses distribution in initial design disk for (a) horizontal and 
(b) vertical orientation 


Fig. 11. (a) Von Mises stresses distribution in initial design disk from centrifugal loads; 
(b) Disk's displacements in axial direction. 
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The results for the stressed state of initial design of the disk from centrifugal load are shown, 
in Fig. 1, Maximum von Mises stress în the disk from centrifugal load is at least 12 times 
greater than in the case of loading from flat of wheel. 


3.4 Shape optimization of cross section of the ellipsoidal disk 


The shape af cross section of the disk is optimized, taking into account only centrifugal load, 
The constructive restrictions allow changing the disk cross section shape only at one side 
and in radial direction at range 150 mm to 300 mm from ellipse center. The section of the 
‘disk at radial distance 0-150 mm has constant thickness b= 10 mm. Three methods are used 
to define the cross section shape (Fig. 12): a) with NURBS knot points, b) with NURBS 
Polygon points and c) with points that are connected with straight lines. Four parameters 
are stated to define the shape. Parameters are varied in the following ranges: 4< X1 £ 10; 4 
© X28105 X35 12,35 X45 for variants a", “e” and 3£ X1s10;05< X2<20;5< 
X3625;25 X425 for “b” The design of experiments is calculated with MSE criterion for 4 
factors and 70 trial points, This design of experiment is also available on the web: 
http://www. mmd.rtu.lv. 


So the 70 strength studies are calculated for each considered method. SW Simulation results 
(volume, maximal von Mises stress, axial displacement of the disk etc.) are entered into 
EDAOpt for approximation and subsequent global search. 


Some optimization and approximation characteristics are shown in Table 4, Results of 
variants "a" and “b” are obtained with second order local polynomial approximation, Third 
order local polynomial approximation is used for variant “c”. Gaussian kernel coefficient a 
was varied for least value of crossvalidation error (4). 


[Approdato] Volume z Maimai von Miss stress 
tk ne [MPa] 
Variant | a w Imm) pe 
Eror Eror 
Auns | o | Predicted | reat | EROF [predicta] reat [F 
8) 1%] 
z as [u [Tamar o] o Sam pe a 
P ua am | oaa [oarran | oms | ao faasa] a77 


E 


seiso [overs | nom | 39998 [2125750] 505 
i ae — [pmo] a A | 


Table 4. Quantitative data of approximation and shape optimization of the ellipsoidal disk 


The obtained metamodels are used for optimization of factors. The ellipsoidal disk volume 
is minimized by taking into account the specified constraints on displacement and stress 
level. The obtained shapes are presented in Fig, 12. As shown in table 4, the best results are 
obtained for variant “b” (Fig. 13), where the volume is lower by 82 % in comparison to 
variant “a” and on 2.6 % -to c". All 3 variants give significant advantage in volume (28.1 — 
33.9%), comparing to variant “d" the initial shape design with constant 10 mm thickness, 
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Fig. 12. Results of optimization of ellipsoidal disk. Shape of cross section is defined by 
(a) NURBS knot points, (b) NURBS polygon points, (c) points that are connected with 
straight lines. 


Fig- 13. Half of 3D model of optimal shape disk and real von Mises stresses distribution in it 


3.5 Summary 
The proposed equipment allows using the standard wheel pair with removable 
measurement equipment as a tensometric wheel pair that considerably reduces material and 
time expenses required for preparing testing. By means of size and shape optimization, the 
total volume of the mounting disk of railway vehicle measurement system is reduced by -64 
‘in comparison with the initial design. The method based on NURBS polygon points gives 
the shape with at least 3% better objective (volume) than other used methods, 
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4. Automotive vehicle gage panel 


In this section designing the mechanical part of an automotive vehicle gage panel is 
discussed. Automotive vehicle gage panels (GP) must meet many requirements - such 
functional characteristics as appropriate stress levels under loads, eigenfrequencies, 
stifinesses, weight, accuracy etc. and last but not least they must have minimal 
environmental pollution during service lifetime. The 3D geometrical models of the gage 
Panel are elaborated using SW. Static and dynamic responses of the gage panel are 
Calculated using SW Simulation and impacts to environment are evaluated using SW 
Sustainability that include such indices as total energy consumed, carbon footprint, air 
acidification and water eutrophication. The stationary and transient behaviors of the gage 
Panel under dynamic excitation as well as stress distribution under static loading are 
investigated. Due to the complexity of the gage panel FEM models, the appropriate 
metamodels are elaborated based on design of experiments, These metamodels are used for 
multiobjective optimization using a global search procedure. Partial objectives are 
aggregated in the complex objective function for optimization purposes. Dynamic behavior 
of the gage panel is then verified by solution of the full FEM models in case of random 
vibrations 


41 Specific requirements, 


‘A constantly pressing problem is the development of safe and environmentally friendly 
engineering objects with high functional properties, attractive style and competitive price. 
We should try to take into account not only precisely measurable functional indices, but also 
such a difficult-to-formalize index as style of GP. 


The Industrial Designer Society of America defines industrial design as the professional 
service of creating and developing concepts and specifications that optimize the function, 
value, and appearance of products and systems for the mutual benefit of both users and 
manufacturer. In fact, industrial designers focus their attention upon the form and user 
interaction of products. There are five critical goals (Ulrich & Eppinger, 2008): 1) Utility: The 
product's human interfaces should be safe, easy to use, and intuitive. Each feature should be 
shaped so that it communicates its function to the user. 2) Appearance: Form, line, 
Proportion, and color are used to integrate the product into a pleasing whole, 3) Easy 
‘maintenance: Product must also be designed to communicate how they are to be maintained 
and repaired. 4) Low costs: Form and features have a large impact on tooling and 
production costs, so these must be considered jointly by the team. 5) Communication: 
Product design should communicate the corporate design philosophy and mission through 
the visual qualities of the products. The practical concept selection methods (Ulrich & 
Eppinger, 2008) vary in their effectiveness and include the following: 1) External decision: 
Concepts are turned over to the customer, cient, or some other external entity for selection. 
2) Product champion: An influential member of the product development team chooses a 
concept based on personal preference. 3) Intuition: The concept is chosen by its feel. Explicit 
criteria or trade-offs are not used, The concept just seems better, 4) Multivoting: Each 
member of the team votes for several concepts. The concept with the most votes is selected, 
5) Pros and cons: The team lists the strengths and weaknesses of each concept and makes a 
choice based upon group opinion. 6) Prototype and test: The organization builds and tests 
prototypes of each concept, making a selection based upon test data. 7) Decision matrices: 
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he team rates each concept against pre specified selection eriteria, which may be weighted. 
The concept selection method is built around the use of decision matrices for evaluating 
each concept with respect to a set of selection criteria. At the same time such formalized 
methods are elaborated as method of imprecision (Zimmermann, 2001) with non- 
compensating aggregation and compensating aggregation as well as fuzzy design method 
with different level interval algorithms, 

The GP styles of different cars significantly differ and should be evaluated in context of the 
specifie vehicle. At the same time style determines an arrangement of particular components 
(distances between gage axes ete: In Fig. 14 we can see initial styles and the 3D model af 
the GP designed for new vehicles (Company Amoplant, 2011). 


Fig 14. Frontal view of GP of the initial styles and 3D geometrical model of G 


4.2 Vehicle GP optimization problem 


Now we can try to use the previously discussed approach for multiobjective shape 
‘optimization ofthe GP. The problem is tated as follows: 


min F(s) =F (a) FF OP © 


subject to g1) £0, j= 1,2. m, and hfs) =0,1=1,3 -€ 
‘where k is the number of objective functions F 

m is the number of inequality constraints; 

eis the number of equality constraints; 


is a vector of n design variables, 


First let us briefly discuss obtaining responses of GP for each particular objective 
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4.3 Strength calculation of GP design 


Generally the strength of the GP design is checked on special vibrostands, The GP is 
subjected to different dynamic loads. Vibrostability and vibration strength of the GP are 
checked on excitations in the frequency domain from 10 to 250 Hz. One of the main 
natural experiments is a test of shock resistance of the GP design under acceleration level 
a = 10g. Such experiments require significant material and time expenses and for 
‘optimization purposes the computer based design check must be used. The 3D 
geometrical model (Fig. 14) of the GP is created using SW and it consists of 18 parts: 6 
deformable bodies and 12 rigid bodies that take into account the inertial characteristics of 
the internal devices. The deformable parts are made from the ABC 2020 plastic, but for the 
internal device bodies are assumed as alloy steel. The initial volume of the GP assembly is 
Do = 764674 cm and mass my = 1.02 kg. The 3D model of the GP assembly is used for FEM 
analysis by SW Simulation to evaluate different responses of the GP. The FE mesh (Fig. 
15) is generated with curvature based mesh (max elements size = 9 mm, min element size 

8 mm, clement size growth ratio = 1.5), that ensures accurate discretization of the 
complex shape bodies of the GP. The FE mesh consists of ~ 210,000 nodes, ~147,000 
elements, ~640,000 DOF, 


In the initial design of the GP von Mises stresses from impact loading are shown on Fig. 15. 
We can see that maximal stresses are concentrated on the brackets cross-section and it 
reaches 4 MPa, Other parts of the GP design are stressed considerably les. This implies that 
the bracket design should be improved. 


‘ig. 15, Meshed 3D model and von Mises stresses distribution in initial design of GP. 


44 Frequency analysis of GP 
Frequency analysis is implemented to find natural frequencies of the GP model and evaluate 
possible resonance in the case of external excitation, The same FE mesh for model as 


considered before is used. The contacts between assembly's parts are defined as bonded, 
The numerical solver FFEPIus of SW is used for calculations 


The obtained results show that the fundamental frequency of the GP is sufficiently high fi = 
17047 Hz. The obtained mode shapes for the GP natural frequencies (f= 201.35 Hz, f= 264 
Hz, f= 331.85 Hz) are shown on Fig. 16. 
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Fig. 16, The mode shapes of the four lower natural frequencies of GP. 


4.5 Sustainability analysis of GP 


SW Sustainability allows getting immediate feedback on the carbon footprint and other 
environmental impacts of the GP throughout its entire lifecycle, including material selection, 
production, transportation, use and end of life (Fig. 17) 
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Fig, 17. Environmental impact of the frame component of GP calculated by SW Sustainability. 
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4.6 Shape optimization of GP 


In this specific situation one of the solutions could be increasing the GP bracket cross-section 
thickness and changing the shape at the mast stressed place. The cross-section shape for 
brackets strengthening is defined by the 3 knot points (Fig. 18 a) of NURBS, Design 
parameters are coordinates of the knot points varied in the following ranges: 3 X1<6; 2< 
X255;0¢ X3 <3. Asa cross-section profile is defined, the 3D-shape is created using the 
Path curve (Fig, 18 b). The same shape for strengthening is created on the second bracket of 
the GP frame component, A maximal von Mises stress in the bracket was minimized with 
constraint on the GP volume (0 <770500 cm3). 


ves 
no a E 
© 


© © 
Fig- 18, Shape definition of bracket: (a) cross-section shape; (b) 3D- shape creation through 
Path curve; (c) shape optimization result 


Von Mises stresses are compared in the design of the obtained shape and the initial shape 
(Fig. 19). There are 6 check points that show von Mises stresses distribution in the most 
stressed bracket cross-section. Volume of the obtained design is» = 770430 cm’. Change of 


w o 


Fig- 19. Von Mises stress distribution in considered cross-section of: (a) initial design of GP 
and (b) optimized design of GP. 
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the assembly volume is insignificant, but maximal von Mises stress level is reduced on -82.4 
Von Mises stress level changes in the cross-section of the GP bracket for intial and 
optimized variants are presented in Fig. 2 


Fig. 20, Von Mises stress distribution in the bracket check paints. 


Fig. 21. PSD of vertical displacement in characteristic points of the GP. 
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Of course, there are possibilities of further optimization of the path curve shape (Fig. 18b) 
and taking into account simultaneously the additional particular objectives such as maximal 
‘Yon Mises stress in material of the GP from dynamic loading in case of harmonic vibration 
excitations; styling of the GP using a-cut method (Zimmermann, 2001), the carbon footprint 
and other environmental impacts of the GP throughout its entire lifecycle calculated by SW 
Sustainability, as well as natural frequencies of the GP. 


For the obtained optimal solution the dynamic behavior of the GP must be verified in case of 
the uniform hase random excitations (see Fig. 21) by analysis of the full FEM model. 


5. Conclusion 


The results of shape optimization of the mechanical components for two different 
measurement systems are presented. The described approach allows obtaining smooth 
shapes that are easy to implement technologically. The jagged forms are excluded from the 
optimization process and there's no need for excessive computational resources. The most 
time-consuming step of the current approach is the FEM analysis of the full model for 
variants defined by design of experiments, the results of which are used for building 
metamodels of appropriate responses, Then the solution of various single objective 
Problems and the implementation of different aggregation strategies for multiobjective 
‘optimization are relatively easy in order to obtain an acceptable final solution. 


It will be interesting to compare the effectiveness of the current approach by using the 
metamodels obtained by kriging and radial basis functions instead of locally weighted 
polynomial approximations 
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